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We study low-lying multipole strength in neutron-halo nuclei. The strength depends only on a few 
low-energy constants: the neutron separation energy, the asymptotic normalization coefficient of the 
bound state wave function, and the scattering length that contains the information on the interaction 
in the continuum. The shape of the transition probability shows a characteristic dependence on few 
scaling parameters and the angular momenta. The total El strength is related to the root-mean- 
square radius of the neutron wave function in the ground state and shows corresponding scaling 
properties. We apply our approach to the El strength distribution of n Be. 

PACS numbers: 21.10.-k, 24.50.+g, 25.20.-x, 25.60.Gc 



Electromagnetic dipole strength in stable nuclei is con- 
centrated in the giant dipole resonance. In recent years 
it was found that there is an appreciable low-lying El 
strength in light neutron-rich nuclei [lj, |2|. This phe- 
nomenon is especially pronounced in neutron-halo nuclei, 
where it can be explained essentially as a single-particle 
effect 0,0. Examples are n Be H© and 19 C Q. This 
low-lying strength was observed mainly in Coulomb dis- 
sociation experiments, for a recent review see [8J. 

The dissociation cross section directly reflects the 
dipole strength distribution or reduced transition prob- 
ability dB(El)/dE. Photodissociation of and radiative 
capture leading to halo nuclei provide equivalent informa- 
tion. The measured strength distributions show simple 
features. They tend to be universal when plotted in the 
appropriate reduced parameters [9j because low-energy 
processes do not depend on certain details of the interac- 
tion. Consequently, effective field theories are nowadays 
used for the description of halo nuclei [10j. In Ref. [llj 
nonrelativistic two-body scattering by a short range po- 
tential was studied using the renormalisation group. It 
was found that the expansion around the nontrivial fixed 
point is equivalent to the effective-range (ER) expansion. 

It is the purpose of this Letter to apply the ER ap- 
proach to one-neutron halo nuclei and find scaling laws 
for the transition strength. This work is similar in spirit 
to effective field theories: the effects of unknown short 
distance behaviour are parametrized by a few low-energy 
constants. It is not our aim to relate these parameters to 
the many-body physics. 

General expressions for electromagnetic strength dis- 
tributions in halo nuclei were derived before (see, e.g., 
[12) and references therein for proton+core systems), 
however, often neglecting the nuclear continuum interac- 
tion. Also expansions for small relative energies were ob- 
tained 13] . The application of ER theory to the study of 
halo nuclei was considered before, e.g., in the description 



of ground state properties |14| , and for radiative capture 
cross sections into s-, p- and d-w&ve bound states taking 
into account the interaction only for s-waves assuming 
a zero-range potential [15j . The effect of the interaction 
in the continuum states on direct neutron capture and 
photodisintegration of 13 C was studied in [16j and |17| . 
While they find a sensitivity on neutron optical model 
parameters for s->p capture, this sensitivity is strongly 
reduced for the cases of p — » s and p — > d capture. 

The prototype of a neutron halo nucleus is the 
deuteron. Radiative capture or photodissociation at low 
energies arc well described by the binding energy and the 
effective range, see e.g. 18]. Quite similar, the low-lying 
El strength in neutron-halo nuclei is determined by two 
parameters if the continuum interaction is neglected: one 
is the binding energy (or neutron separation energy S n ) 
which determines the overall shape of the dB(El)/dE 
distribution, the other one is a normalization constant. 
In the analogous case of 19 C, the binding energy could 
be determined directly from the shape of the observed 
strength distribution [7j . Effects of the neutron-core in- 
teraction are parametrized efficiently by the ER approx- 
imation. 

The cross section for photodissociation of a nucleus a 
into a neutron n and a core c is given by 
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The reduced transition probability for electromagnetic 
multipolarity nX is denoted by dB(nX)/dE. The photon 
energy is given by E 1 = E + S n with S„ — h 2 q 2 /(2fi) 
where q is the inverse decay length of the bound state 
wave function and /i is the reduced mass. Using de- 
tailed balance, the corresponding radiative capture cross 
section can be obtained from this expression. With the 



equivalent photon number, cross section JI} determines 
the (first order) Coulomb dissociation cross section |g. 

The strength distribution dB(EX,li — > lf)/dE for a 
certain transition depends on the corresponding matrix 
element that contains a radial integral with the wave 
functions of the initial and final states with orbital angu- 
lar momenta U and If, respectively. (We consider spinlcss 
neutrons and cores, a generalization is obvious.) At low 
energies E — h fc 2 /(2/i), the main contribution to the 
radial integral arises from radii larger than the radius 
R of the nucleus. This will be true for all possible val- 
ues of the angular momenta involved. In the figures of 
Refs. [lalUl it can very well be seen that the radial inte- 
grals are dominated by the outside region. Away from a 
resonance, a case which we always assume here, the con- 
tinuum wave function is small inside the nuclear radius. 
In a hard sphere model it is exactly zero. 

For the neutron+core case the dimensionless reduced 
radial integral is well approximated by 



q. Clearly, for small separation energies of the neutron a 
large transition strength can be expected. 

The reduced radial integrals (J2J) have the general form 



#W 



7exp(— 7) / 7 \ '/ 



( 7 2 + K 2)A+ 



r© 



(5) 



H { + )lf (A) cos(k + S tf ) + n { ~ )li (A) sin(/t + 6 tf ) 



with k — kR where the rational functions 1Z 
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be found for all relevant values of U and If by means of 
the recursion relations 
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(h\ ') functions, respectively, that describe the behaviour 
of the radial wave functions beyond the range of the nu- 
clear potential. The phase shift Si, contains the infor- 
mation on the interaction in the final continuum state. 
Introducing the dimensionless shape functions 
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the reduced transition probability is given by 
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with the effective charge number Z^' = Z c [m n /(m n + 
m c )] A and the asymptotic normalization constant (ANC) 
Ci i of the bound state wave function. For neutron+core 
systems the dipole transitions are dominant since El 
transitions are suppressed by the small effective charge. 
The ANC of the actual many-body wave function as it 
appears in J3J consists of the ANC of a neutron single- 
particle wave function and a spectroscopic factor which 
accounts for the many-body aspects. The single-particle 
ANC is determined by the normalization of the wave 
function. It depends on q with the scaling behaviour 
\C U \ 2 oc q for l t = and \C k \ 2 oc qj 21 ^ 1 for h > as 
seen, e.g., in a square- well model. Here, we have intro- 
duced the dimensionless scaling parameter 7 = qR. The 
factor |CzJ /q 2X+3 in the reduced transition probability 
(ff|) depends only on properties of the ground state. It di- 
rectly shows the scaling with the characteristic parameter 
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with spherical Bessel (ji), Neumann {yi), and Hankel starting from 
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The radial integral involving the regular scattering wave 
function shows a r A -';+'/+ 1 behavior for small values of 
r, thus it can be extended to R — 0. On the other hand, 
the radial integral involving the irregular function shows 
a r x - l i- l f behaviour. Thus for A = 1 we can extend the 
integral to zero only for li = 0, If = 1 and li = 1, If = 
and have convergence. In the other cases, one has to cut 
off the integral for a finite R value. 

For low energies we can use the effective-range approxi- 
mation k 2l+1 cot(5i) = —1/ai + rik 2 /2 + ... for the phase 
shifts Si with the scattering length a/ and the effective 
range n in order to take the interaction in the contin- 
uum state into account. (Note that ai and n have the 
dimension of a length only for 1 = 0.) We restrict our- 
selves to the lowest order term, i.e. 
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with the dimensionless reduced scattering length q de- 
fined by c 2 +1 = ai/R 2l+1 and the ratio x = ft/7 = k/q = 
\jEjS n . We assume that this is a reasonable approxima- 
tion of the phase shift. We further assume that the scat- 
tering length has a "natural" value with ci being O(l). 
For example, for s-wave scattering from a hard sphere 
with radius R we have Co = 1. In exceptional cases, this 
assumption is not fulfilled, e.g. if there are resonances in 
the low-lying continuum or subthreshold states. For an 



s-wave halo nucleus bound by a zero-range force the scat- 
tering length is given by clq = 1/q, i.e. cq = 1/7. In such 
a case the scattering length in the s-wave is unnaturally 
large; it diverges in the halo limit 7^0. 

We can expand the analytical results for S t f (X) in 
terms of the small parameter 7 and obtain for the most 
important cases for dipole transitions 
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The term in front of the square parenthesis in Eq. IjlOjl 
is the well-known result without continuum interaction, 
see, e.g., [5lll7|. In general there will be a modification 
of the shape function. It is getting more and more im- 
portant for larger binding energies of the halo nucleus, 
i.e. larger 7. Transitions with If = l t — 1 are more af- 
fected by this final state interaction than transitions with 
If = li + 1. In the case s — > p the first correction to the 
shape function for 7 — ► appears only in the 7 3 term 
proportional to c\ . There are no corrections of lower or- 
der in 7 independent of the final state interaction. This 
explains the remarkable uniformity of the characteristic 
shape for the low-lying dB(El)/dE strength in halo nu- 
clei like 11 Be and 19 C. For the p — > s transition there 
is a greater sensitivity to the final-state interaction. It 
appears already as a correction linear in 7 and cq. This 
is in accord with the results of Ref. [ljj . 

The total transition strength can be obtained by inte- 
grating @ over the relative energy E. It is given by 
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For dipole transitions and no final-state interaction (i.e. 
q = corresponding to a plane wave in the final state) 
we find 
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A finite value of 7 leads to a reduction of the total 
strength in the continuum as compared to the extreme 
halo limit 7 = 0. With final-state interaction, the func- 
tions T l f (A) also depend on the reduced scattering length 
ci, however rather lengthy expressions are obtained. 



TABLE I: Dimensionless quantities c[, see Eq. ©, and spec- 
troscopic factors C 2 S from the fit to experimental data [g] of 
11 Be Coulomb dissociation. 



model 


3/2 

Cl 




c\ /2 C 2 S 


effective-range 
Woods-Saxon 


-0.41(86, 
-0.46(70, - 


-20) 
-14)" 


2.77(13, -14) 0.704(15) 
1.87(3)" 0.696(15) 



"At energy E = 0.4 MeV in the peak of the strength distribution. 

The value of B(E1, k) = J2i ■ B(E1, h -> l f ) for all El 
transitions from a state with orbital angular momentum 
li is directly related by the non energy-weighted sum rule 



B{El,k) = 



Z (1) e 



2_3 

4tt 



— (r 2 )i- 



(17) 



to the root-mean-square radius (r 2 )i t of the same bound 
state, see, e.g., 0. The rms radius scales with 7 for dif- 
ferent values of the bound-state orbital angular momen- 
tum k as (r 2 ) ex R 2 h" 2 , (r 2 )i oc R 2 /j, and (r 2 ) h ex R 2 
for li > 2 1^. On the other hand, the energy- weighted 
sum rule (Thomas-Reiche-Kuhn sum rule) gives a con- 
stant independent of U. From the above we see clearly 
that the low-lying strength is most pronounced for s-wave 
bound states and to a lesser degree for p-wave bound 
states in halo nuclei. Considering the q dependence of 
B(El,li — > If) in (fni)! the scaling of {r 2 )i i with 7 is ex- 
actly reproduced. With no final-state interaction the full 
transition strength as predicted by the sum rule is found 
in the continuum. However, if there are bound states 
in the particular final state the continuum contribution 
to the total strength is reduced and the shape will be 
distorted. This will correspond to a non-zero value of q. 
Let us now apply the general theory to n Be, which 
shows a very pronounced single particle halo structure 
with a neutron bound by 504 keV in the i ground state 
that can be considered as a 2s\/2 state. Assuming a 
radius of R = 2.78 fm the neutron separation energy 
corresponds to a small value 7 = 0.4132 < 1 typical for 
a halo nucleus. Unlike the case of the deuteron, where 
there is no excited bound state, there is a \ -bound state 
at 320 keV excitation energy. We can regard it to a 
good approximation as a lpi/2 single-particle state. We 
perform a y 2 fit of our analytical ER approach (not using 
the expansion for small 7) to the experimental Coulomb 
dissociation data 13 by adjusting the reduced scattering 
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lengths c/ and c/ in the two p-wave channels (see 
Table HJ) and the ground state ANC Co . The quantity 
C]/ is definitely unnaturally large as expected from the 
existence of the weakly bound | state. In Figure ^ our 
model fit (solid line) is compared to the experimental 
strength distribution 6] , including corrections due to the 
detector response. 

The obtained ANC C = 0.724(8) fm" 1 / 2 can be trans- 
lated into a spectroscopic factor C 2 S by comparing to 
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effective-range model 
Woods-Saxon model 
transition to 1/2 bound state 




FIG. 1: Reduced transition probability as a function of the 
excitation energy E* — E + S n compared to experimental 
data extracted from the Coulomb breakup of 11 Be 6]. 



the ANC of a normalized single-particle wave function 
generated from a Woods-Saxon (WS) potential with a 
depth adjusted to the correct neutron separation en- 
ergy. Assuming a radius of R = 2.78 fm as above 
and a diffuseness parameter of a = 0.65 fm we obtain 
C 2 S — 0.704(15) which compares well with the spec- 
troscopic factors from experiment and theory as given in 
Ref. [6( . The model ground state wave function has a rms 
radius of 7.172 fm that corresponds to a total strength 
of B(E1) = 1.645 e 2 fm 2 according to (fT?}l . Applying 
the spectroscopic factor we obtain a total strength of 
1.159(36) e 2 fm 2 . The experimental continuum strength 
6] exhausts about 78(5) % of the sum rule. The ground 
state transition probability to the first excited state is ex- 
perimentally known to be B(El) = 0.116(12) e 2 fm 2 
corresponding to 10(1) % of the sum rule. This bound 
state transition strength is indicated in Figure ^ as a 
dotted line assuming an arbitrary FWHM of 0.160 MeV. 
The remaining 12(6) % can be attributed to thep 3 / 2 state 
that is occupied. Comparing the experimental transition 
strength to results of a plane-wave calculation would un- 
derestimate the extracted spectroscopic factor. 

Our ER approach can be compared to an exact calcula- 
tion of the matrix elements with wave functions obtained 
from WS potentials. We assume the same shape of the 
potential with R and a as given above. Repeating the 
X 2 fit by adjusting the depth of the potential in the two 
p waves independently we find C 2 S — 0.696(15) close to 
the result of our ER approach. The exact calculation 
(dashed line in Figure^ almost coincides with our ER 
approximation. We can also extract the dimensionless 
quantities c{ according to Eq. 10. Their modulus de- 

1/2 

creases with relative energy (much more strongly for c/ 
than for cy ). Since the c{ in our ER approach are con- 



stant they can be considered as an average of the actual 
values over the peak of the strength distribution. Tak- 
ing this into account a reasonable agreement of the two 
approaches is found again (Table QJ. 

In conclusion we find that low-lying electric multipolc 
strength in one-neutron halo nuclei can be described ef- 
fectively by a few low-energy constants: the binding en- 
ergy, the scattering length and the asymptotic normal- 
ization constant. We apply our ER approach to the halo 
nucleus n Be and find a remarkable agreement with an 
exact calculation. The exhaustion of the strength due to 
the bound state affects the extraction of the ANC and 
of the spectroscopic factor. In our approach there is no 
need to determine parameters of optical potentials from 
elastic scattering. 

The present theoretical approach will be applied to 
other halo nuclei. It will be extended to higher multipo- 
larities and to proton halo nuclei in a forthcoming pub- 
lication J20J. An effective-range theory of two-neutron 
halo nuclei like n Li would be very interesting but is cer- 
tainly much more complex [4|. An early work related to 
this problem is |21|. 

We are grateful to T. Aumann and H. Emling for useful 
discussions. 
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